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T Mies thesis research was conducted as part of a 
continuing project at the Naval Postgraduate School (NPS), 
which is investigating the effects of the addition of )isiimag 
particulates to solid rocket propellant to enhance the 
performance characteristics. In order to correctly Wage. 
the effects of the particles, an accurate measurement of the 
sizes and quantities of the burning solid particles must be 
made with relation to their axial position within the port 
of the propellant grain. 

Previous research has attempted this data collection 
through a number of different methods. [Ref. 1] This 
thesis concerns itself only with the use of holograms of the 
bUrENine propellant! in. a cross flow section of the combustion 
chamber. The methods of obtaining these holograms at’ the 
NPS are discussed in detail in References 1 and 2. 

Once the holograms are obtained, the problem is 00a 
analyze the data they contain. Considerable effort has been 
taken in this area to utilize a Quantimet /20 image 
processing system and a PDP-11/04 computer to conduct 
automated data retrieval on the reconstructed holograms. 
Tis research has been the subject of two previous’ theses 
[Ref. ivand 37 which desctibe the process. This method of 


image reconstruction and data formation has experienced 


IG 


Sermerakemajyor dlificultie’s. 

The Quantimet 720 has been the largest obstacle to this 
method. Memnas demonstrated very low reliability, with the 
system constantly unusable due to intermittent faults. Data 
outputs have varied significantly between sequential runs on 
the same hologram images when the system has operated. 

Fniewadditcion — the Quantimet /20 requires that a clear 
distinction exist between the grey levels of the background 
and the object features in the input image. The system is 
unable to carry out any image processing to separate these 
two levels itself. The holograms used inherently contain 
light speckle which varies over the complete range of grey 
levels in the hologram. The Quantimet has been unable _ to 
distinguish between this speckle and the real features. The 
Besult is that the background and the features’ have 
overlapped and true feature counts have not been able to be 
ebtained. (Ret. 1 and 4 | 

The Quantimet results required initiation of research 
into alternate methods of data retrieval for the holographic 
images. Pnompat. eu la tr; digital processing was felt to be 
required to separate the image features from the background 
of the hologram before sizing and counting the features. A 
set of ImageAction software was purchased to see if this was 
feasible. As Em 1s software Et Seasonal IBM/AT computer 
augmented with a TV image acquisition and display plug-in 


board, the concept of carrying out the complete image 


ital 


enhancement and feature collection on that system appeared 
og real: In this manner the unreliability of the Quantimet 
system could be completely overcome. 

Hence, rather than continuing with the previous (you 
with the Quantimety = 207 this thesis deals with a new 
approach to the problem. The initial requirement was _ to 
evaluate the ImageAction software to determine if it was 
able to reduce the background sufficiently to permit gees 
extraction from the holograms. With this capability sien 
hand, the next step was to write a software algorithm which 
was able to retrieve all the required feature data from the 
enhanced image. 

This thesis is divided into four further  ¢hajpieeuae 
Chapter iT desc 1 bes eine ImageAction software, its 
capabilities and the routines of main applicability tosiimes 
problem. A discussion of the format that the software uses 
to store images is also carried out, as it presented an 
obstacle to the further processing of the enhanced image: 
Chapter III contains a description of the program which was 
written to perform the feature data extraction one 
enhanced image. Chapter IV describes the results obtained 
by using the completed IBM/AT processing method to extract 
data from a standard calibtatiwon Vrecie we Chap temas 


contains the concluding remarks. 
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lig Ba HOmOGhAtEMAGE PROCESS ING 


The ImageAction software acquired for this project 
required both hardware and software to be installed on the 
IBM/AT. The hardware consisted of a PC Vision Frame Grabber 
board which plugs into the IBM/AT computer. This board takes 
an analog input video signal and digitizes it; assigning 
each pixel an integer value. (mismameanmtal data is then 
stored in the frame memory on the board and is available for 
PRocessitte or transfer via the computer bus. The stored 
image can also be viewed, as the display logic on the board 
transforms the digital pixels back into analog form before 
Semaine it 6a video monitor. The PC Vision board and the 
ImageAction software were both made by Imaging Technology 
mon po raeed MO Carry out digital image processing. 

The software package consists of a comprehensive, menu- 
riven set of routines which can carry out image processing, 
image graphics and limited image analysis. Once installed 
and loaded, ImageAction appears on both the computer screen 
and on the image display monitor. The monitor screen is 
used to display the image being processed and _ also the 
current software menu in use. The computer screen shows the 
software status, current prompts for both keyboard and mouse 


Satries, and a help block if required. 


slags 


From the main menu one can access via the mouse, al l Yorn 


the submenus to carry out the required po rd@eqanunce F i eilgee 


2.1 shows a complete breakdown of the - ImageAction menus. 


(Reto 








f il fle 
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Figure 2.1 The ImageAction Menus 
[Ref. Sj jURiecwrer 2-2) 


The ImageAction software permits 


the processing of the 


entire MOnitor screen or it allows the subdivision of the 


screen into four quadrants. It is also possible to designate 


any rectangular portion of the screen as an Area of Interest 
3 


14 


CAGT)., anni to move this area into any quadrant for easier 
processing. Reno ener very useful algorithm permits the entire 
screen to be compressed into one quadrant of the screen. 
iis process is called a Quad aaa 

Pema eneay compe Lronrm this function ES important as it 
allows a full screen image to be processed and later have 
the data extracted in one fourth of the time. However, it 
must be noted that in order to do this, tne algorithm only 
copies the even-numbered pixels on the even numbered lines. 
This implies that some of the smaller image features may be 
lost and that the resolution of larger objects may _ be 
affected. 

Once an image is processed, graphics can be added to 
further enhance the image or to record specific information 
@n it. The graphics Sa EERE i e's imerudce titling and script 
text as well as the drawing of limited geometric shapes. 

With this general overview of the ImageAction software 
immeis MOw possible to discuss its detailed application to 


the problem of hologram image processing. 


fee ACQUIRING THE HOLOGRAM IMAGE 

The ImageAction software requires a stored video image 
mp Ut. This means that the hologram must be reconstructed 
and then viewed by a video camera system, whose output 


either is fed directly as an RS1/0 standard video signal to 


5 


the PC Vision board or is stored 6n a Video medneaeee that it 
can later be played back into the POC Vision bia 

Figure 2.2 shows the hologram reconstruction  t¢€c hitagme 
A more detailed discussion on the actual mechanics Of Me 
process is available in Reference 1. Ini tiger, a VHs 
Camcorder was interfaced with the microscope. Later the VHS 
camera was replaced by a number of Newvicon cameras with 
increased light sensi ae A discussion of the Camels 


used is included as Appendix A. 
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MICROSCOPE 


REF EWEN 
ANGLE 


OBJECTIVE 
LENS 


RECONSTRUCTING 
KRYPTON. [ON 
LASER BEAM 


Figure 2.2 Hologram Reconstituc ureon 
[Ref .. Pyare cea 
The microscope can be set Gkor a variety of 
magnifications by rotating a front turret ring to aie cee 
desired set of optics. Images can be acquired using 


microscope objectives with magnifications of one, two, foams 
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or ten. Meiteoamer Levels Of Magnification are used, the 
fopemeotmireld becomes smaller and so the focusing of the 
image becomes more difficult. Also as the magnification is 
miereased, the field of view becomes narrower and so eee 
light is available. Hence, more sensitive cameras are 
required. 

In addition to the hardware shown in Figure 2.2 for the 
reconstruction and-viewing of the hologram images, 26.15 
possible to place a spinning mylar disk at the location of 
the reconstructed real image. This disk serves to reduce 
the speckle in the reconstructed image but it also reduces 
mee amount Omi echteawhichseenters, the optics. This loss’ of 
intensity was originally a problem that was later overcome 
by using a more sensitive camera. 

The video output of the cameras could be directly input 
into the PC Vision board as previously stated. However, due 
to the separate physical location of the reconstruction 
laser, the images were stored on a standard Beta video 
recorder and played back into the board. 

Once the image is stored on video tape, there are three 
meenods for entering it into the board. The first two, Snap 
and Grab, are similar in that they take one frame of the 
recorded image and display it for processing. The Grab 
technique is preferable as it permits the played-back image 
to be viewed until the desired frame is shown and then, by 


use of the mouse, that frame is “grabbed”. 
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The third method of acquiring an image in the )compuee 
is called Average. This technique takes a specified numoen 
Of frames from the video @oinpweeeanma averages thei. 
intensities, producing a final, single time averaged {eau 
As the speckle background varies with time, this techna@ie 
of acquiring the image was of some interest and was 
investigated in detaiel. 

With the hologram image now in the computer, the Wenge 
step was to attempt to enhance the image by reducing the 
background. The goal was to completely remove the speckle 
background from the image while not losing any of the real 


feature data. 


ee IMAGE ENHANCEMENT 

The PC Vision board assigns a grey value between zero 
(which is black) and two hundred and fifty-five (whichees 
white) to each pixel of the image. In an ideal image,the 


particles to be analyzed would be of one grey level, 


preferably dark, and the background would be of a separate 
grey level, usually lighter. Figure 2.3 shows a histogram 
of the grey levels of an ideal image. The grey scale is the 


horizontal axis and runs from zero on the left to age 
hundred an fifty-five on the right. The vertical axis is a 
measure of the number of pixels of each grey value present 


in the image. 


i 


impeach icemme ene background of an image is never all of 
Gomeo rey levebedud Melthner are the pixels of the particles. 
Even in a very clear white light reconstruction of an image, 


there are some bands of grey levels which exist for both. 





Puree 2e5o Histogram of an Ideal Image 


Figure 2.4 shows a histogram of a clear white light 
reconstructed image which has both particles and background 
present. (No speckle is present in the background.) The 
particle grey levels are to the left and are clearly 
separate from the band of background pixels to the right. 
Thus, if a threshold level was set between the two bands, 
all of the particle pixels could be set to one value, while 
the background pixels could be dhe be another. The ideal 


case could then be realized and the image could be analyzed. 


Ne, 








Figure 2.4 Histogram of ayWhitemla wnt nace 


The hologram image background, however, is varied in 
ithe s 1Cy « It is generally more intense (lighter) than the 
desired particle features, but the speckle has salt and 
pepper intensities that can range over all of the grey 
levels present in the image. Hence dark speckle can appear 
on what should be a white background giving the appearance 
of more particles and lighter speckle can appear in the 
particle images causing them to appear to have holes. 

Figure 2.5 shows a histogram of a hologram image (ieee 
has particles preseme. The effect of the speckle is that 
there is no clear band of grey values for either the 
particles or the backerounde This makes the assignment of a 


thiresho id value impractical as the speckle exists 
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throughout. Hence the image cannot be converted back into 
the ideal case and the sizing analysis of the image cannot 


be done. 





lg ahem tblye = es) Histogram of a Hologram Image 


What is therefore required is a system to process out 
Mies speckle ain the image in order to establish a clearer 
separation between the grey levels of the particles and 
those of the background. The Quantimet system has no image 
processing capability and so cannot remove this speckle. 


iimes itapility has greatly limited its use in this research. 


[Ref. 1] 
lee Filtering 
The ImageAction software has many methods ot 
filtering and image processing. Of interest to the hologram 


Taga 


image problem are the algorithms tor convo tmemom These 
routines compare each pixel to the surrounding pixe Vo 
change the center pixel's value based on a weighted kernel. 
In this manner very small unwanted speckle can be convolved 
or f£iltere@—Mouwpeot the backoaound. It must be noted however 
that this also tends to blur or smear the a@ctual partie 
shapes resulting in some loss of edge resolution. 

Two of the ImageAction routines are particu lae 
useful for reducing the presence of speckle. These are the 
Lowpass routine and the Blur routine. The kernel designs for 
these two algorithms are shown in Figure 2.6. [Tt can .eige 
seen that the Lowpass routine involves ay three-by— ime 
convolutrone while the Blur involves a five-by-five 


CONnVOLUt1oOn. 
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Figure 2.6 Convoltition Kerne bey aiues: 
(Ref. 5] 

If the speckle is spatially very small, the Lowpass 
routine can increase (ie., brighten) the speckle grey level 
without blurring the edges of the particle features 
significantly. However, if the speckle is larger (as is the 


case in the higher microscope  magmiticasemonmmom the 6B iim 
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roOUbpune muctebe Used to completely remove the speckle. This 
Serr clOme recs mecmime slenificant smearing of the particle 
shapes. 

While neither of these routines provides a wide 
ape bathom Ol wl éCatires jana background to @locate a threshold, 
use of these routines with a calibrated particle array can 
effectively remove all of the speckle and provide a useful 
threshold value as described in Section 4.B. 

The ImageAction software also permits user-designed 
convolution kernels and several other weightings have been 
tried. It has been found that the results of the Blur and 
the Lowpass algorithms are sufficient for the problem _ to 
this date. 

Zs Dhoee sine bd ins 

Once the background speckle has been sufficiently 
removed by eee | nce the entire image can be thresholded to 
completely remove all of the background. Tet Som ly 
required that the filtered background be increased to just 
one integer grey level more than the grey level of the 
brightest particle image. (Further separation causes a 
Eurther Pmereascrtn resolution. } Hence while some speckle 
may be still very apparent to the eye in the filtered image, 
it will be totally removed by thresholding. 

The thresholding algorithm assigns a lower and an 
upper threshold value. All pixels which fall inclusively 


between these two values are give a value of zero and 
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therefore become blacks All grey values outside of the 
range are given a pixel value of two hundred and £ itty =e 
and are therefore made white. 

With the image in this form it is now possib commu 
extract the feature data information of the paneeemeumen 
without any interference from the background. ImageAction 
has no sizing capability to carry out this extrac 
unlike the Quantimet. Also the ImageAction software does 
not permit another program to bDesrun concurrently. Thags 
implies that the processed image must be saved as a data 
file, so that another program can access it to Carry OWteeee 


feature data extraction. 


C..- SAVING THE ENHANGED IMAGE 
The ImageAction software permits the processed images to 
be saved to- disk. The image is stored in a string file that 
has a sixty-four byte header, followed by up to thirtyeeue 
bytes of comment and then the image. Assuming that the image 
has been compressed into one quadrant by the Quad Squish 
routine to speed up the processing, the image is 65585 
bytes long. 
_ it must be noted that this image file is a string “foe 
that can be over 65,600 bytes long. Fortran does not easily 
read string sequences-7 and the length of this striae 


prohibits the use of any standard word processor to reformat 


al ear 
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MpeeGtrder sLOemake the Stored image file easily readable 
Boa POrtianmpmoezrzam, 4a short Basic program was written to 
convert the file into a column file with a fixed 14 format. 


The details of this program are included in Apperidix B. 


Presse ULOMATIC PROCESSING 

The ImageAction software also permits repetitive image 
processing techniques to be programmed to run automatically. 
Mimics facility is done by the Script menu functions. Once 
the script function has been initialized with a file name, 


all subsequent steps are recorded until the script function 


mee turned off. Then, by simply entering the subroutine's 
name while in the correct script sub-menu, the recorded 
steps will proceed automatically. This allows the hologram 


images to be rapidly processed and stored with operator 
interaction only when required. 

With the image processed and stored in a Fortran 
readable format, Dietow is possible to proceed to extract 


the particle feature data from the image. 
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pis PARTICLE DATA COMPIRAT ION 


It is important to recail that the? intention of eae 
research was to produce a system which is capable 
measuring the quantity and size of solid particulates in the 
combustion chamber Ofmwawrockercamotcoce These particulates 
are expected to range in size from 1 micron to 200 microns 
in diameter. [Ref. 2] They are also expected to be nearly 
spherical in shape. 

The system must therefore be capable of counting the 
number of particles present in a given area of the 
combustion chamber. Further, it must be able to Size themmuan 
terms of area or diameter. It must be able to distingumem 
which particles are round, since this capability can be used 
to exclude particles which are of irregular shape (such “as 
fuel binder and ‘celratft), and also to remove multiple 
overlapping particles which would upward prejudice the size 
data. With these requirements as the basis, a Forttang 
program was written to carry out the data extraction "tam 


the processed image. 


Aw THE PROGRAM FORMAT 
The program was written in Fortran for several reasons. 
Birs &. hans language could easily accommodate all of the 


required functions and routines. Second, if this program was 
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PomEDCUcedeimwconjUunction With the Quantimet /20, the 
Muentimete?20.allso uses Fortran as its programming language. 

Third, the ImageAction software has been introduced as a 
stand-alone driver package called ITEX/PC. Eneet basso rm, 
each of the functions in the menu-driven format can be 
called separately as subroutines. This .ariver sathkeae is 
available in the Fortran, Pascal and © languages. Hence, 
these subroutines ,could be easily imbedded into the data 
Sx ELaActlomo m= program to create a single, complete algorithm 
for both the image enhancement. and data tabulation, when the 
driver package was purchased. (This work is part of a 
follow-on thesis effort) 

The program was written in a sectioned form, with each 
major step being written as a stand-alone routine. In this 
format it was felt that if any of the ImageAction driver 
routines were able to carry out one of the functions in a 
more economical manner, then it could be easily inserted. 
Also, this permitted each of the steps to be tested for 
correctness and for economy of time. This becomes important 
minmen Processing up to 65,600 pixels of information. 

The program was written to analyze a stored and 
reformed image file of one quadrant in size. This implies 
that the image to be processed must either have been only 
one quadrant large to begin with, or the full image must be 
Quad Squished prior to storing. Additionally, the image must 


have been stored with no comment included. If a comment is 
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inc Luded, the program will read the comment as part of Segme 
image and all of the data will be incorrect. 

It was recognized early in the image processing portion 
of the thesis that each image was different in terms of 
background intensity, particle distribution and usable image 
area. This was in some cases due to the camera used, but 
was also a function of the smoke in the background and _ the 
alignment of the recording and reconstruction lasers. Hence, 
it was realized that the data extraction algorithm would 
have to be able to accommodate these variations. Tike 
program was therefore written ina fully interactive format, 
with the user being prompted for all required input 
parameters or being given the option of using some 
predefined values. 

Also included in the program is a predefined test image 
which can be used at any time to ensure that the algorithm 
iS  -workKing “correctly. This is particularly useful (wime@ 
changes are made to the routine. The test image and _ the 
correct outputs are shown in the comment statements of the 
DrOgram Prior to each sSub=-fLunerion:. 

A complete copy of the program is available from _ the 
Thesis Advisor. Due to its length, it is not included Tenge 
B. PROGRAM BREAKDOWN 

As stated, the program is broken into sections. Each 
section will be briefly described, with the program 


Limitations outlined. 
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dhe Image Input 


Micweeiis (eae rtom Fregulred 1s tO input the image to 
be processed. The user is offered the Gmorece of using the 
test image or a data file image. If the data file image is 
selected, the user is then prompted to define the actual 
area of the input image which is to be processed. 

The variable frame is defined in terms of the X and 
mmecoOrdinates of the upper left and lower right corners of 
the desired area. If the entire image is to be processed 
then the user selects 0,0 and 256,256 for the two corners of 
the variable frame. 

The data file is then read in. The program 
Pifoma cical ly strips off the first 64 data entries as these 
are the ImageAction header bytes. It then stores the image 
mmeo a 296 by 256 array. It must be restated that the image 
must have been stored without a comment block, for if it has 
a comment, that comment will immediately follow the 
ImageAction header. Thus, the comment will be read in as 
part of the image and the image array will be destroyed. 

ikea the variable frame area is less than the full 
frame, the image is read into the top left corner of the 
array. All subsequent processing will only be done on the 
used area of the array. 

The user is then given the option of viewing the 
image to ensure that it has been read in correctly. The 


image is shown on the computer screen, with each pixel being 
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the three digit integer value of its grey scale Teo) eee 
the top left corner of a large image can be viewed, however, 
this is normally enoteh@te ee if the image appears to be 
COLLeC.. If the image appears incorrect the uSer Can /Sitauee 
the routine over again. 

a Image Thresholding 


The image is now thresholded. The user can define 


the upper and lower bounds of the threshold region eae 


Lf ibea values are usually determined by interactive 
thresholding while in MImageAction (before storing the 
filtered image). The array is reduced to zeros fone 


pixels outside of the threshold region and to onesie 


pixels inside the region. Once again the user is given the 
option to view the thresholded result. If the result asi 
as desired, the user cameestar t Gere rc. It should be noted 


that the thresholding permanently alters the pixel array and 
so if the thresholding is not correct the image must Deueee 
entered. (The program automatically reprompts the user for 


this, Tequirement.) 


De Feature [dentification and Labeling 
The thresholded image 1s now converted to a 
feature-labeled image. The program does this by working 
through the image a row at a time. starting in the “ippen 


left corner, each new pixel in turn is evaluated. {1f i tae 
zero it is skipped. If 1t ds. ap ones the pixel to its left 


is evaluated. If the left pixel is non-zero then the new 
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Dimcilme is piven the same value as its left partner. If the 
left pixel is zero then the pixel immediately above the new 
Peaelice testedewith similar results. Should both the pixel 
to the left and the pixel above be Zero, then the feature 
count register is advanced one and the new pixel is assigned 
Enact number. This sequence is continued until all of tthe 
rows have been labelled in this manner. 

Once all of the rows are done then the program goes 
through each row again, checking each pixel with the - one 
epove iit. Tienes cuUrrentupaxe!l isezero it is skipped. ag 
the pixel above it is zero, or of the same value, it is left 
Math its current value. If the pixel above has a different 
value then all of the pixels with the value of the current 
pixel are changed to the same value as the pixel above. In 
this manner features which have a He We in them, or a section 
which is partially detached, are sure to be assigned the 
same feature number. 

Once this row-to-row check is done, then all of the 
features are renumbered from one to the total number. The 
final result is that the array now contains a fully feature- 
labelled image, which can be processed for specific feature 
data. The user is once again given the option of viewing 
this image before continuing, and the total number of 


features is displayed. 
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4. Feature Statistics 

The feature labeled image is now used to find ieee 
specific feature data which is required for the ho lgmeeae 
image analysis. This includes each feature’s X chords 
the horizontal diameter), the Y chord (or the Vey tmeem 
diameter), and the feature's total area. Each of these 
measurements is done as a pixel count. 

The X chord is found by counting all | GfRsae 
continuous pixels of the same feature value in each _ row. 
Should a feature have a break in it, additional sums foupmeme 
same feature are maintained for the row until the end of the 
row 1s reached. The largest sum is then used as that “moms 
chord length for that teattire: Then that row's chord values 
are compared to previous row's chord values for that feature 
number. The largest X ee ee is kept, and the program 
eontinues to the mext Tow. 

The feature areas are computed at the same time as 
the X chord lengths by simply summing the total number Saige 
pixels encountered in each row for each feature number. 

The Y chord lengths are computed in the same manner 
as the xX chord lengths with the exception that the values 
are calculated column by column instead of row by row. 

These measurements are all that is required to 
assess the hologram image for its particle content, as_ the 
other wiactors., such as circularity, can be found Erom eens 


Although all of these measurements are in terms of pixels, 
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memwoumge=npe cCasy to convert them to actual size measurements 
DC cyOwrercrence pixel size was available for the 
magnification of use in the image. 
3) - Feature Data Display 

Piwopde: Lomcacny OUlLeamalysis of the feature data, 
it must be displayed in some usable form. The program was 
eee en toeworter the data in two formats, depending on its 
intended use. The user is prompted for the desired form of 
output. Poem a ches PlaveeOresElve Lt eature data, as either a 
Gable=~or a histogram, the user is given the opportunity to 
select additional data outputs, either in the alternate form 
Sein a modified version of the current form. 

a. Piseogtan Cutput 

If the histogram form is selected, the user is 
advised of the total number of features and the maximum X 
euordgd length of all of the features. Prompts are then 
issued to define the lower X chord limit of the histogram, 
the upper X chord limit of the histogram and the number of 
Memeo that the Hhistogram is to be broken into. 

Next the user is asked if non-round features 
are to be excluded. If they are, a test flag is set and the 
program proceeds to form a histogram file. 

The feature data is then sorted into the bins 
With each feature being counted into the bin of the correct 
emord size. Each feature is tested for roundness, eg the 


moet tae is set, by finding the magnitude of its X chord 
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less Ltsoe Y chlo rd. If this magnitude is greater thane 
then the feature is excluded as non-round. The value Of ia 
was found by experimentally checking a number of hologram 
images and finding that over 95 % of the features me t sue 
criteria. In fact, 85 % had a magnitude less thane. 

The bin values are then converted 1 ite 
percentage values by dividing each by the total number of 
features. The final result is then displayed as a data file 
and the user is advised of the total number of features and 
the number of features which were excluded due to non- 
roundness. 

The user is then asked if a hard copyaiiame 
requested. If one is, a copy of this histogram data 1S iitima 
copied to” (disk, where it can be later printed out or up- 
loaded to the NPS mainframe computer where it can be 
converted to graphical form by a plotting Tomeunee 

Be Fabdular .Ouep wt 

If the tabular output is selecrcdr then tie 
feature data is printed out by feature number, displaying 
the area, the X chord length and the Y¥ chord lenge The 
user is again asked if a hard copy is required and, if )Sonea 


copy of the table is written to disk. 


Cr PROGRAM LIMITATIONS 
The program was written for use on the IBM/AT and 


because of this, suffers from two main limi tae temse F i posite 
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the memory space available on the IBM/AT and the program 
sizes are such that images with a maximum of only two 
Mimd@eedmand titty, (250) particle features can be processed. 
The program will run with images that have more than this 
Number, however, the feature specific data is lost. 

iiemseconmdsmajor; limitation is the time that it takes to 
run this program for large image sizes. Lindt quadrant 
is to be analyzed, _the program takes approximately two and 
mpe=-nhalt hours to run from start to finish. This is due to 
the number of operations and the speed of the IBM/AT. 

BOLE OLetmese problems could be overcome by running the 
program on the NPS mainframe computer. However this would 
require that the stored image be uploaded to the main 
computer and this in itself requires considerable time as it 
Ho a file of over 65,600 bytes. MeLUGener difficulty would 
be encountered if ImageAction drivers were to be included in 
the program, as the hardware and software for their use does 
not exist on the mainframe computer. 

CPewelsebeltm tliat "a better alternative is to fully 
investigate the use of the ITEX/PC driver subroutines on the 
BEM / AT. It is possible that with their use the program may 
be substantially increased in speed. Eta is  .algsio ENO 2 hit 
that additional memory will be freed as the pixel array may 
moe be Pees. since the PC Vision board can be used to 


store the image during processing. 
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D. SUMMARY 


The feature data extraction program has been used on many 


processed images and has been found to operate well, given 
its) two wmelimieationce The use of the LITEX/PC  “Gisaieauee 
SUbTOUti nes May greatly enhance its performance. iL oe 


additional measurements or calculations are requaned® they 
can be easily added into the routine due to its section 
format. . 

The program was used to produce the preliminary results 


which follow in Chapter sa. 
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ees Payot wna RY RESULTS 





This chapter contains the results obtained by using the 
MmageAction software and the Fortran data retrieval program 
Pomeowenminedmam the earlier chapters of this thesis. The test 


images for all of the results in this chapter were produced 


meom a calibrated standard reticle. 
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Pustren see Calabration Standard Reticle 
[Ret.06, Pigure 1 | 
imtcemmneELcCheCNwCONSiIsStS Of Cwo sections as seen in Figure 
eel ice Cimeathhar sipOGtdon) Of the reticle consists of 
approximately 10,000 opaque circular features of twenty- 
@onmee We esizesenancinge frometive (5) to ninety=three (93) 
fiecrons. These features are randomly photo-etched into the 


Pee (oO) millimeter circular area. The array to the right 
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of the eight millimeter circular test pattern contadmsueeee 
twenty-three standard sizes of features arranged in rows of 
five features as shown in Figure 4.2 

The calibrated array was used to investigate the effects 
of the image processing on the size of the proces 
particles. The manufacturer's data for the particle jeu 
distribution of the circular reticle portion of ° the stanuauen 
calibrated image was used to evaluate the results obtained 


from the IBM/AT system. 
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Figure 4.2 Calibration Array Sizes (Niezonw 
| Riewee 201 

This calibrated standard reticle was available in two 
iO ins. It was purchased etched onto a glass plate which 
could be used to carry out white light analysis Of jem 
system. It was later recorded onto a laser hologram in the 
Same manner that the holograms were produced of the 
combustion products in the rocket motor. It should be noweg 
that thrs test reticle hologram was recorded without the 


smoke present that occurs in the combustion holograms. Thus, 


Sie, 


a discussion of the effects of smoke on hologram analysis 


Satepot be carried out. 


a. THe beeeto OF IMAGE AVERAGING 

i Wes Deak! ec womope mites uexhibited by lasem reconstructed 
hologram images have been previously discussed. As this 
speckle pattern is a random property, which can be 
considered time-varying due COmVvVibretromawand the laser 
coherence, miomwas @ielwmtinat a lareeuportionwof this speckle 
could be removed by time-averaging the reconstructed image. 

In chapter two, two alternative methods for time 
averaging of the reconstructed image were presented. These 
two methods cco the use of the Average routine of the 
ImageAction software, Seen em ine lWSton Gf va Spinning mylar 
Mesk in front of the microscope optics. The relative merits 
of each of these techniques was investigated and the results 


mon low: 


ies ImageAction Averaging 
The ImageAction Average routine takes successive 
image frames and averages their grey levels. In this manner 
the random time-varying features in an image will be washed 
Sut, White the fixed features will not. Usame this routine, 
image lengths of one, two, bawe:, ten, fifteen and twenty 
frames were averaged. The resulting images were then 


fmertcered and thresholded to see the effects. 


oes 








An Unaveraged Hologram Image 


Figure 4.3 





Twenty Frame ImageAction Average Hologram Image 


Figure 4.4% 
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The results were as expected. The averaging 
greatly reduced ene speckle, even to the unaided eye prior 
Pomtiltering and thresholding. An example of an unaveraged 
hologram image is shown in Piers To eetneexample of a 
twenty frame averaged hologram image is shown in Figure 4.4. 

It should be noted that both of these images have 
not been filtered or thresholded. A comparison of the two 
figures shows that the time averaging has removed or reduced 
the lighter speckle out of the solid features. It has also 
reduced the amount of darker speckle from the lighter 
background. 

It was also found, as expected, that by averaging a 
greater number of frames that the speckle was more reduced. 
This reduction occurred up until ten frames were averaged. 
The results for averaging een frames and twenty frames 
were found to be the same for that of averaging ten frames. 
aus, Pte wes £LOund  that™an application of the Average 
routine on ten frames was optimum. 

or. iy Pere Averaging 

When the spinning mylar disk was inserted into the 
reconstruction equipment several interesting results were 
noted. ole gs ae the amount of light incident on the camera 
was greatly reduced and a much more sensitive camera was 
required. However, Hitcum proved snot to be a problem as 
cameras were available which could easily make up for this 


Woss (as discussed in Appendix A). 
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It was also observed that the speckle was greatly 
reduced. te is due to the fact that the mylar @iskigiee 
rotated at a rate greater than the 1/30 second intéeeracue 
time of the video camera. The mylar disk, therefore, changes 
the phase of the incident laser light over time, and emme 
changes the speckle pattern seen by the camera. Hence, the 
speckle values become washed out while the constant particle 
images do not. 

Once again this was evident even to the unaided eye. 
However, when the image was filtered and thresholded in the 
same manner as the image using the Average routine, it was 
found that the two final images were identical. Hence the 
spinning mylar disk carried out the same function as_ the 
Average routine applied over ten frames or more. 

The image oat by the spinning mylar disk was 
available immediately from the reconstruction apparatuse 
While the time-averaged ImageAction image required a time- 
consuming averaging process. Therefore, the spinning mylar 
disk was felt to be the better solution to achieve a time- 
averaged image. 

Or Mylar Plus I[mageAction Avetaring 

in interest of completeness, an ImageAction 
Average image was made from one, two, five and ten frames of 
a mylar reconstructed image. Even after filtering and 
thresholding no improvement of the speckle reduction was 


found, even over the ten frame averaged image. Thus, the 
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averaged image produced by the spinning mylar was found to 


EPemCOpeinal wimeporn speekle Teduction and in processing time. 


Breeetbre ERP PBCTS sO IMAGE PROCESSING 

The hologram images, even after time averaging, Ste) 
exhibited a substantial speckle content. An example of this 
can be seen in Figure 4.4 and oe oD he Ure 4. shows 
the same time-averaged image as in Figure 4.4, only now 
thresholding has ‘been done to necneant eenes speckle “that 


remains in the picture that may be undetected by the eye. 





Figure 4.5 Time-Averaged Hologram Image after Thresholding 


Therefore the image must be further processed in order 
to remove this speckle prior to the feature data being 


extracted. WemerwitsSewethiscs speckie will be Seotnted as 
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Features. It is this remaining speckle that was ueeeeuee 
partial cause of the Quantimet system's inability to 
correctly perform feature extraction on the recons triemaee 
image. [Ref. 1] 

The two further processing operations which are regimes 
are filtering and €hresholdimne. These two operations were 
described in Chapter Two. The effects and limitations mie 
these functions will now be covered. 

Las Piltenime 

The two predefined ImageAction types of convo /ijGaam 
are Blur and Lowpass as described in Chapter two. Lf time 
size of the remaining speckle after time averaging is large, 
then @the Blur function must be Weed: By ‘large’, ities 
meant that the speckle appears to be of the same size as the 
smallest of the known feature sizes. In cases where the 
speckle is not large, the Lowpass function is preferable as 
it causes less smearing of the edges of the feature 
particles. In general, it was found that the Lowpass 
function could be used for microscope Magnifications am 
times “twomand Less. 

When the correct filtering routine was applied, all 
of the remaining background speckle could be completely 
suppressed. However, it should be restated that the smaller 
particle features are also being suppressed by t hase 
filtering as are the edges of the larger features. Also, ie 


two particle features are close together, their edges may be 
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PutimrucdmmtOete ther to form one solid particle. thas” last 
peopmemel smpareieularly evident at dower magnifications. 
Me Thresholding 
With the speckle suppressed by ae Tatas bee et ie mete 1S 
now possible to produce a thresholded image Eheaee Camara 1 ts 
wl yesknown particle features and no Hackzground. An example 


of this is shown in Figure 4.6 





Peaitem4.O wmeawerazsed, Pidtered and Thresholded Image 


ou BPmocessine Lamitations 
Several bitte smeCane be quickly noted Erom Figure 
ee 0. Tt can be seen that this image does not show the 
Sopebnre Calibrated image as is shown in Figure 4.2 . The 


particles of diameters 5, (Ls on, a 17, and 21 have been 
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lost. However the background which is "evident ea Fi pee 4.4 
and Figure 4.5 has been completely suppressed. 

It was also evident that the size of the particles 
which were present in the image, could be varied in sizemeg 
adjusting the level of the threshord- The image shown in 
Figure 4.6 shows the calibrated features at the largest 
threshold value before background speckle begins to appear. 
However the meaning of the pixel size of these features is 
clearly dependent on this threshold level. In a mixed 
target such as a combustion hologram, the question arises as 
to where to set the threshold level, as it will be hamwda tie 
tell what is speckle and what are small particles. 

The most obvious solution to this problem is to include 
the calibrated reticle on the back of the combustion chamber 
so tha€ it becomes part of ecaé@¢h hologram Wamacer The 
threshold level for the hologram can be determined from this 
image and then applied to all subsequent images from that 
hologram. This solution is dependent, however, on a CONS tame 
background average grey level, which may not be the case in 
a smoke-filled combustion chamber. Other techniques shoul 
be investigated to remove the smoke or other spatitag 
variations of transmissivity if they become a problem. Some 
alternate techniques include bleaching the hologram [Ref. 4] 
or background grey-scale averaging techniques. [Ref. 7a 
This determination of the chaeenone setting, however, lis 


clearly a limitation on the system. 
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Mitlietemitemeis Mecessary to consider why the six 
lest DpanrtEtclc f@atiires were missing Im the processed 
hologram image of Figure 4.6. Looking at the images before 
Gncesmoldine (Figures 4.4 and 4.5) it can begseen that these 
six features do not readily appear in either of these two 
Figures. It may be possible that these small particle images 
were UWntresotved in the holographic recording of the 
calibrated reticle and they were not removed by the 
Bilteriwng andmenresholding. 

Figure 4.7 shows a white light image of the 
ealibraced array before it was recorded on a hologram. (at 
is felt that the physical damage to this array that removed 
some of the deposition occurred after the hologram was 
made.) The smallest particle features can Deen Seen Tuten is 
Peeure “at the arrows). Figure 4.8 shows this same image 
after it has been filtered and thresholded. The three 
smallest features of Figure 4.7/7 have been removed by the 
Peocessimg operation. 

Hence, several reasons for ‘the loss of the six 
smallest features lve the hologram reconstruction are 
possible. Firsts some of the features may not have been 
mero ve@dmimethe recording of the hologram of the calibration 
array. This could be due to a variety of reasons. It may be 
that the holographic process is unable to resolve particles 
eof this size, however, Shrsebss mot felt to be true in this 


case. 
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SS 





Figure 4./ White Light Calibration Array @mace 
Showing the Three Smallest Features (Arrows) 





Figure 4.8 Filtered and Thresholded Calibration Array 
Showing the Three Smallest Features Removed 
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AlMeernatively, tiemiremocram reconstruction process 
may cause particles of these sizes to be unresolved. This 
could be due to many reasons such as the intensity of the 
moconatructnonmemase Don ame Oorrect focusing of the optics, 
MPieorrect alignment of the optics, insensitivity of the 
mecording cameras, etc. All of these areas need further 
investigation. 

iW sticmmmtrtre =bteht reconstruction, only the top 
MuiGeemmparticeles became unresolved after the filtering and 
m@iresholding (Figures 4.7 and 4.8). However, the next three 


Sezes Of particles were still clearly visible, S@motlteri ng 


should not affect these particle images. Therefore, in the 
hologram case the remainder of the particles should have 
survived the processing, however, an additional three 


particle images were lost due to unexplained causes. 
EUS: ree that the small array particles 
became unresolved due to a combination of the above causes. 
Regardless of the causes, the reduced resolution of the 
processed hologram image is a limitation on the system. lise 
this limitation cannot be overcome, then the effect of the 
the inability of the system to reconstruct and process 
particles of these smaller sizes must be assessed with 
respect to the goals of this data retrieval system. 
meeey ORT FICATION OF THE IBM/AT SYSTEM 
Patne wehe MMinmitations of the system into account, the 


entire IBM/AT hologram data retrieval system was tested on 
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two calibrated téstiiimacese The two images used were a 
2X magnification of the white bien reconstructed 
calibrattionm retwtele, and a 2X magnification of ) the are 


reconstructed. hologram of the calibrarewon sree moe 
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Figure 4.9 Manufacturer s Size Distr bamerog 
[Ret "Gl 
The manufacturer's particle size distribution forme 
calibration reticle is shownsingecure 6 Jee Each of (ttre 
particle sizes is plotted as a bar to show i1ts re alee 
density. In order for this data to be compared to the two 
reconstructed image results, it is required that the 
diameter sizes be converted from lengths in. Mienonse fo 


lengths in pixels. The details of this conversion asmveue as 


a0 


the truncation of the manufacturer's data to exclude the 
unresolvable smaller particles are contained in Appendix C. 
iD Winste bieht Calabration Reticle Results 

Die manufacturer s particle size Aiguripmedene have 
been replotted for the white light reconstruction case. The 
resulting histogram is shown in Figure 4.10. As previously 
stated, details on how this modified plot was made are 
available gin Appendix Ce It is assumed that this data does 


not include any overlapping particles. 
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Pigune 4 miamumacturer s Size Distribution 
NGiiwecmhtome sReconstruction 2X Magnification 


The white light reconstructed image was taken using 
the spinning mylar disk as the time averaging device. Dit 


was then lowpass filtered, and the resulting image was 


eee 


thresholded to the Level required to just remove Gime 
background speckle from the calibrapionma ae The image 
was then stored in accordance with the procedures outed 
in Chapter 2 and Appendix B. 

The Fortran program of Chapter 3 was then ~cUree 
extract the particle “dataz The tabular output andre 
histogram data output are included as Appendix D. 

The “plotted histogram of the results is showniieea 
Pere meas.’ < It can be seen that the plot compares well 


with that of the manufacturer's data in Figure saneee 
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Figure 4.11 White Light Experimental Size Distribution 


The general shape of the distribution of "Su 


manufacturer's histogram and the experimental data histogram 


DZ 


Ne ———————————— 


are the same. However, the experimental data indicates that 
more of the smaller particles have survived the processing 
Bone prwedictea., i Lhws could™be due to a variation in the 
level of the background intensities between the area of the 
array and the area of the reticle pattern. This would 
result in the threshold level being set too low and some of 
the speckle features may have been included. 

Also a very small number of particles appear to be 
larger than expected. Pits ls spr opaply due to overlapping 
particles which passed the roundness test. However, the 
mun DeLwor patrtacles in this group is small enough that they 
can be neglected. 

hus, Lat appears that for the white Pieht 
teconstruction case, this method of processing and data 
recovery works very well with the exception of the smallest 
particle images. 

oe Laser Hologram Calibration Reticle Results 

The manufacturer's particle size distributions have 
been replotted for the laser hologram reconstruction case as 
described in Appendix C. Since the unresolved particles are 
Maceer for this case, the curve is different from the white 
Bret case. The resulting histogram is shown in Figure 4.12. 

The laser hologram reconstructed image of the 
calibrated reticle was taken using the spinning mylar disk 
as the time averaging device. It was then processed and 


thresholded in the same manner as the hologram calibrated 


Sys: 


array. The thresholded image was then stored and the data 
file converted in the same manner as the white lien cme 
The Fortran program was run to “extTraicas the 
PartLrclhe weadata. The tabular output and the histogramemaemee 
output are included as Appendix E . It can be seen that 
there were more particles present in the hologram imageuseuem 


the white Light image. The plotted histogram is Shown 
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Figure 4.12 Manufacturer's Sizem@Drstr rumen 
Laser Hologram Reconstruction 2X Magenitieatran 


Again it can be seen that the plot compares well 
with that of the converted and truncated manufacturer Ss eqeusee 


in its general size distri pw@eren- However, a larger number 
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of small particles are present in the experimental curve. As 
Piemewidth of the one pixel diameter range is now 13.1 
microns instead OF feo ni Gr ons:, ba ies felt that these 
patiti¢c les “are still speckle which has Survived the 
processing. 

Ptwecan also be Seen that a ftew larger particles are 
miesomonce Again present and it is still felt that these are 


maerto Overlapping, particles. 
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De sUMMARY 


The results shown in this chapter have been obtained by 
using the IBM/AT system. They. have shown that the system 


can produce usable data results from hologram images, with 


8 


some stated limitacionmse: The loss of small particle images 
when using holograms requires further research to determine 
which particles are lost at each step of the process. La 
this manner corrective action can be attempted to overcome 
all “or some Of thewloccece The effect of those (Sean 
particle sizes which cannot be recovered will have to be 
determined. 

The effects of uneven background illumination caused by 
the presence of smoke in the rocket chamber also remains to 
be investigated. It is clear that if this smoke is presume 
then it will have to be removed by some process. If ite 
small enough it may well be removed by the filtering and 


thresholding. This remains to be seen. 
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V5 CONCEUSLONs 


The IBM/AT system for the processing of laser 
reconstructed hologram images has been discussed and _ the 
Poawliminatyeeresults from this system Nave been analyzed in 
Miapter LV.=) fhe combination of the ImageAction hardware and 
Beimeware with  thet= Fortran data extraction program has 
Deoamueced results which compare favorably with those expected 
feom calibrated sources. Therefore, it is felt that the 
IBM/AT system should undergo further research. 

The system has three main limitations on its use which 
need further study before it can be considered truly 
effective. Pie iostearmpOtrtant of wtnese limitations is the 
Mess of the small particle features Po the hologram 
images. Piremmactual == cause Of the loss of all of these 
Particles is not fully evident. 

From the comparison of the white light reconstruction 
example and the laser hologram reconstruction of the same 
Bai bprated array it can be seen that the filtering and 
processing of the image may cause the particles of less than 
mine (9) microns to be lost. Further work may show that 
meet erent Convolution kernels may be able to retain even 
these small features. However, the comparison of these two 


maces also shows that an additional three particles (with 


DH 


- 


sizes up to twenty-one (21) microns) have been lost ameeme 
hologram image. 

The fact that these particles do not even appear tome. 
eye in the unprocessed averaged image, Suggests that they 
may be lost prior to the precessimee Further investigation 
is required to determine if they are lost in the hologram 
recording, in the hologram reconstruction. in the image 
recording or in the-image processing. 

The IBM/AT system can be used to test other techniques 
to improve the first three of the potential loss areas. ie 
may be found that the hologram photographic images do not 
have these size particles inathem- However, once the best 
results have been achieved, the final minimum sS1i1Ze00iom 
particle: YEhat can survive the entire process will have to be 
assessed against the aim of this yeeas 

Assuming that the first limitation of the system GsSeiteue 
to be acceptable or can be improved upon, the second main 
area of concern is the length of time that it takes _ to 
completely analyze an image. At present, once the image is 
stored in ) video. form, it can take up to three hours see 
process one full quadrant of an image. This includes the 
time to process the image, save it to a data got Lee convert 
the data file, run the data extraction program and themepssem 
the resulting data into graphical form. This amount of time 


can definitely be reduced. 
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The subroutines that make up the ImageAction software 
Vimeo vatlabple as andependent drivers in the ITEX/PC 
software set. It is felt that the use of these routines 
should be fully investigated to try to speed up the entire 
IBM/AT system. Dieetnesemrolmemnes can be Called directly 
maome the Portran data extraction program, it is felt that 
much more than half of the time currently required can be 
eliminated. 

Pmatily ene tiinrd matin limrtation of this system is the 
small area that can be analyzed at one time. This is due to 
the lack of memory space on the IBM/AT, as images of only 
250 features can be processed. Oneernacdinmewilt is felt “that 
this problem can be largely overcome by using the 
ImageAction subroutines im the ITTEX/P@ software “directly 
maom the Fortran program. In this manner the monitor can be 
used as the pixel array memory and Mette the 256 by 256 
array in the IBM/AT will not be required. This should free 
up more than enough memory to process an image with as many 
features as are present in a 1X magnification image. 

In summary, this thesis has presented an alternate 
method of performing analysis on hologram images. The 
meMyAL computer is the basis for a digital system that can 
‘carry out image processing as well as data extraction. The 
ImageAction software routines used are discussed. The 
Fortran program written to carry out the data extraction has 


been explained. In the preliminary results the system 


ah 


demonstrated that it has 


the ability to perform the requiem 


functions with some Specific. limi tareone The causes iam 


these limitations are described and possible solutions have 


been offered. 
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BP PENDEX A 


EOWPLLGH tT eCAMERAS 

Lie wetoLocmamemreconstruction process required that the 
image be viewed through a microscope and then be captured as 
a video image. A video camera was therefore required. As 
it was required that the image be stored on tape and brought 
Boma separate location for processing, it was felt that the 
optimum camera would be a standard video camera/recorder 
Ceamcorder) . Phuswethe  tirst camera used was a VHS Movie 
Camera. Detailed specifications on this camera are in 
Kemeotre A.l . 

The VHS Movie Camera provided several advantages. A 
separate recording machine was not required and the camera 
had a baleen background sensing Cae Gulatis which 
automatically adjusted the aperture for varying illumination 
momdi tions . 

Iter however, also presented several disadvantages. 
Parst, the main lens was not removable. This caused 
problems in that the camera optics did not mate well with 
Meme rLoscOpe Camera port. The result was that not all of 
the camera viewing area was illuminated and a black area 
around the outer perimeter of the image was formed. As 
well, due to the number of lenses in series, a darker spot 


was observed in the center of the image. 


Out. 


FEATURES /CAMERA 


Panasonic VHS 


PanasonwTe 


Panasonic 


Movie Camera WV-1410 WV-1460 
ENS 6:1 Power Zoom, None, None, 
F1.4 (C9=stanea Standard Standard 
Not Removable Fol teen Fit time 
PICK-UP 
Te Ee NEWVICON VIDICON NEWVICON 
é 
HOR EZON EAL 525 eNes 6.00 LINES 600 LINES 
RESOLUTION 
MINIMUM 10° Bux -5) Lux O05 Lite 
ELLUMINATLION 


Paecure: Aso Camera Standard Features 


The major disadvantage of the VHS Movie Camera was 
discovered when it was used to view the laser-reconstructed 
holograms. While the camera had worked well for the white 
light reconstructed images, it was quickly discovered that 
the camera did not have the light level sensitivity required 
to record the hologram reconstructions. 

The VHS Movie Camera was only able to be used to record 
images at magnifications of 1X, and then only if (Gime 
spinning mylar disk was removed. It was therefore decided 
that a camera with a lower light sensitivity was required. 

The next camera to be used was the Panasonic WV 1410. 


This camera is for low Light level security surveillance. 


It did not have an internal recording system, and so an 
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external BETA recorder was used. The WV 1410 did not come 
aaneoelcons we but 1t had @ Standard fittine that matched the 
camera adapter of the microscope. His. oncerinstalled, it 
DPaoduced ss images withowt either the outer dark halo or the 
inner central darker spot. 

While the WV 1410 did not have an automatic aperture 
system, the image could be viewed directly by use of a 
ioe ornesena Seehe Light Level of the image was adjusted by 
Varying the power of the reconstruction laser. At the lower 
magnifications of 1X and 2X this camera could be _  over- 
illuminated if the mylar disk was removed. However, it was 
found that at higher magnifications, and at lower 
magnifications with the spinning mylar disk, that the camera 
produced images that were too dark to process. 

Hence the WV 1460 camera was tried. This camera had all 
of the advantages of the WV 1410 camera, and had sufficient 
sensitivity to be overilluminated at up to 10X 
Teter ication with the spinning mylar disk in place. By 
adjusting the laser power it could record images over the 
complete range of magnifications, with or without the mylar 
disk. 

Thus, the WV 1460 was felt to be the optimum camera for 
this system and was used to record all of the hologram 
images used in this thesis. The white light images used in 


the thesis were taken with the Panasonic VHS Movie camera. 
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APPENDIX B 


CONVERSION OF THE STOREDRItivGr Fiver 

The ImageAction software permits an image that has been 
processed to be stored for later use. The program takesmeme 
image and stores it as a string file which contaitnmsuead 
ImageAction header, a variable length comment block and an 
image. bloek. In the image block each pixel is stored as an 
ASCII character with a decimal value from 0 to 25 55qiRaa 
accordance with ts erey ever 

The Fortran data extraction program is designed to  @tame 
this file and read the image portion into a two dimensional 
array, which it can then process. The Fortran prog Games 
limited, due to the 512K memory capacity of the IBM/AT and 
the array is therefore specified to be 256 by 256. This 
corresponds to one quadrant of the ImageAction screen, and 
the software permits a single quadrant to be stored. 

Difficulty was experienced, however, when the data file 
was attempted to be read. Fortran does not permit easy 
access Co string files and, in particular, excessively ome 
string files. If the image is stored without a comment, the 
stored image file is 65,600 characters long. It is ¢c lease 
therefore, that the data file had to be processed before it 


could be read into the Fortran program. 
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Conrvenmtronal ss wotds processors also cannot handle string 
files of the length of this data file. itera ed, C4001, the 
ImageAction header contains some undefined ASCII characters 
and so it was unreadable to them. Using the M/DOS' Debug 
routine it was found that the header was 64 characters long. 
While the Debug routine could have been used to convert the 
image data file, _ this would have been a long manual 
operation, and so this was felt to be unacceptable. 

The final desired form of the image data file was to 
itave st as a fixed format column of data which could then 
easily be read by the Fortran program. Hence, it was decided 
that a program would have to be written in some other 
programming language which could “comwmert.) the file. The 
resulting program was written in Basic and is shown in 
maeure B.1i. 

The program takes the ImageAction stored file, named 
Porn Ge ss and reads a string of data from it. Eteeshou ld bie 
mored thatethe longest string in Basic is 255 characters. 
mayen ASCII character is then taken one at a time and 
converted to its decimal integer equivalent. The integer is 
then formatted into I4 format by adding the correct number 
@it preceding zeros, Ameen itursestored back into the 
Sueput f£ile DATA. The program continues until it has done 
mek 65,600 characters. 

The program experiences difficulty if it encounters a 


character at the end of one of the input strings which is an 
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ASCII control character. The ASCII control character= eh 
integer pixel grey values in the image that range from zero 
(0) to twenty six (26). Dhesescharacten-. force the program 
to abort as they will not permit the string section | )tcgseeee 
closed. This problem can be easily overcome by preprocessing 


the stored image to remove these control characters. 


10 Di EEST G25) ;sVARIABLE DIMENSION 
20 OPEN"O" 51, "A2DA TAS -OUTPUT Rigs 
30 OPEN LI” 72, Ac DA TAs G UE NPU Tease 


40 FOR J = 1 TO 257 

50 LIST1$ = INPUT$(255,#2) :GET STRING OF Dag 
60 FOR I = 1 TO 255 

70°. X$ = MIDS(LIST1I$ em ;SELECT ONE CHARACTER 
80 P = ASC(X$) -GET INTEGER VAMUE 

90 TF P < 1@—cone. a0 

100 IF »P<9100 Goro deo 


110 PRINT#1,P -STORE VALUE OF P 
120 Gono on90 

130 Y$ = SPACE$(2) ;ADD FORMAT SPACES 
140 PRINT#1,Y$;P >STORE VALUE OF P 
150.) Geto 196 - 

160 Y$ = SPACE$(1) -ADD FORMAT SPACE 
170 PRINT#1,Y$;P _;STORE VALUE OF P 
180. TGOmo. 190 

190 NEXT I -DO EACH CHARACTER 
200 NEXT J -DO NEXT STRING 
210. S@iKeS act “CLOSE OUTPUT mite 
22 Ome Ge se 2 -CLOSE INPUT FILE 


Figure B.1 ImageAction Data File Conversion Program 


Therefore, the steps required to store and convert the 
image data file are shown in Figure B.2 . These steps are 
carried out after the image has been processed and _ the 
desired results have been achieved. The image must be 
contained in only one quadrant. Should a full screen image 


be desired for data extraction, the ImageAction Quad Squish 
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Power nentweaempe used Used to compress the image before these 


steps are taken. 


Si i beh Chir 225 
Sob jae: BZ MAP LUT TO AOI 
SEP |S ADD LUT. 30 
a0 Pe 4 WAaEeeuU Tl TO,ACL 
STN 5 es, SAVE AOI AS A:DATA.IMG 
SIS oe. RUN PROGRAM TR1.BAS 
Spa eT RUN FORTRAN PROGRAM 
Figure oe 2 Image File Storing and Conversion Steps 


The ImageAction High Clip routine takes all of the 
image pixels above the stated value and makes them into the 
stated value. The memory values, however, are not changed 
until a Map Lut to AOI operation is done. The top 30 values 
must be removed before the entire image is increased by 30 
because the ImageAction software works as a el loneored Lut 
system (i.e., when you add one to 255 it becomes 0). 

This file conversion process has a large disadvantage in 
that it takes fifteen minutes to run the Basic conversion 
program. It then takes a further eighteen minutes to read 
the converted file into the Fortran data extraction program. 
Memets Celt that use of the ImageAction Driver routines as a 
Dae te of) the Fortran program could possibly eliminate the 
need for many of these steps. This alternative requires 


PuUrther research. 
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APPEND Lie 
CONVERSION OF THE MANUFACTURER Sp 
OF THE CALIBRATION Rat aie 

The calibration reticle was produced by Laser Electro- 
Optics Ltd. The reticle portion of the image ¢onsis tSoeee 
average of 10,490.7 dots which have been etched onto a glass 
backing. For use to testing the hologram imaging equipment, 
a hologram of the calibrated array and reticle was also 
made. The hologram was made in the same manner as_ the 
holograms from the combustion chamber. 

The calibrated reticle consists of twenty-three (23) 
different sized circular particles which are deposi ted@ywagay 
known densities onto the glass in an eight (8.) millimeter 
Gace he. Beside the reticle is an array which contains each 
of the twenty-three (23) sized particles in rows of five. 

Figure C.1 gives a complete breakdown of the sizes of 
the particles and the number of each of them present in the 
peticce we; This data was used to plot the histogram shown in 
Figure 4.9 of the preliminary results. 

The sizes of the particles for the calibrated array and 
reticle are measured in microns. However the sizes of ‘the 
particles measured by the IBM/AT system are naturally 
measured by a number of pixels. A method of converting the 


micron diameters to pixels is required. 
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PART ITCES PEAMETER NUDE ER. 


NO. (Microns) PRESENT (Avg) 
1 5 ele 2904.7 
2 oc 7 eT 
3 8.89 892.9 
4 11.73 1168.6 
5 eee 7 WOO <7 
6 21.46 642.9 
7 7. 00 Hn 
8 26.15 Giles 
9 30.90 397.6 

10 34.40 28,0 

iL 37.12 Be a1 

ee 40.68 F379 

13 43.19 2a 7 

14 47.40 Lore 

iS Sioned 3 106.7 

16 527.32 114.3 

17 56.34 98.3 

18 60.98 85.4 

19 ‘Sy eae | 64.6 

20 73.94 De, 

2A 80.77 iy ae 

22 Be 4.2 

23 92.88 h0 


TOTAL 10490./7 


Rae te GC a Calibrated Reticle Data 
[Ref. 6] 


The size of a reconstructed particle image is a function 
Memmene Macmifteation of the reconstructing optics, the type 
of monitor used, item type of filterimeealcornithm and the 
level of the image threshold. (These is assuming that the 
Mackeround intensity iis constant throughout the image. 
Therefore, any comparison of sizes between images must _ be 


done under these conditions and limitations.) 
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Clearly a white light reconstructed image and a hologram 
image will differ in several if not all, of these aspecue 
and, therefore, each must be considered separately. The 
white light reconstruction will be Consider, cae. 4 

In this case, an image of the calibrated array was taken 
at a 2X magnification and processed. The method of 
filtering and the level of thresholding was recorded so that 
it could be reproduced for the reticle image later. A 
photograph of the white light array is shown in Figure (gee 
The image was then stored and the Fortran program was run to 
extract the feature data. The tabular output has been 
reorganized for comparison and is shown in Figure C.2. 

Using the X width data from Figure C.2 tor alien the 
particles which were present and which were obviously not 
malformed, an average length per pixel was calculated. The 
diameter in microns for each of the particles was taken from 
the manufacturer's data in Figure C.1. The diameter pem 
pixel is shown in Figure C.2 for each feature, as is the 
average. (It should be noted that the data for particle 4% 
was not used in forming the average length per pixel as the 
data had too large a variance to be considered. 

This average diameter per pixel value was then used as 
the conversion factor to convert the manufacturer's particle 
distribu tioms to pixel diameters instead ory micron 
diameters. This method is felt to be a relatively accurate 


Conversion. 
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ils) 4 S, 6220.6 
i 4 5 ais 
17 5 D 6.88 
ZS 2 6 14) 
Ze , 7 8.14 
26 6 6 ea, 
oS, j 8 Oy 7 
4 3 y 8 Us ile, 

Z2 6 A NOTE (1) 
44 8 8 7.04 
igs 9 10 Cec 
74 9 sid 7.4/7 
82 10 ys Tg 

100 10 ile NOTE (1) 
125 ez 14 (EP | 
154 14 es, 6.64 

Roo tG teen 7. So eMC RONG / PIXEL 


Particles were not used as their diameters 
appeared questionable in Figure 4.8 due to damage. 


Particle not included in total average as it 
has too large a variance. 


Figure C.2 Calibrated Array Feature Data 
White Light Reconstruction 
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Using the conversion factor calculated as shown, @ tae 
for the white light reconstruction was made. The Table shows 
the width range in pixels VerTSeS tile Widtherance oy microns 
and the total number of particles which fall into that temas 
range. The number of particles in each range was found by 
summing 


the appropriate totals from the manufacturer 5s =aasee 


in Figure €C.1. The results are shown in the £14785 € tgeeeee 


co. luinins “Of: “Taple €. i: 


TABLE C.1 CONVERTED MANUFACTURER’ S DAT 
WHITE LIGHT REGONSTRUGimon 
DEAME LER DIAMETER TOTAL NO. OF TRUNCATEDSRGS 
(PIXELS) (MLR ONS) PARTICLES OF PARTIGEES 
O - 1 OR OLS = io SOS.6 «4 O=n0 
T° fee a2 1o3) = Laer, 2306 el ces Lise. 
ne) 1G) =e 1Gi0e 16303 
3 - & 22 0 C=) 29a oS 937 
G - § ZO de Se 3 07D 20 60785. .46 
Sia 16 SO = Gar 0 77 Aa ie sla (| 
6 - 7 44.0 - 51.3 2/79 2a 2 7 Fomet 
Tous 51.3 G= wee / Ze 2 2 
8 - 9 ao eae S| 8) en 9, Se 85.4 
oe AO 66 5010=, 74.20 90s 9028 
gE 0 ees ei 1G. 0 = 6 ee 12.4 1 2eo 
Lie 1 Z S1.3 = 985..0 eeea 4 wz 
eZ ere? eS SSO) ro oe 1.0 130 
A histogram of the manufacturer's size 
distribution based on pixel diameters has been plot tedaeag® 
iS shown. ln tae ce C2 oe 
The large distribution of particles in the zero to one 
pixel area of the histogram is only true, however, if all of 


the particles which are on the reticle can be viewed. This 
momeno tthe case Eor the IBM/AT system. The smaller 
papeeeles are mot reconstructed and so they should not be 


mcluded ine the histogram. 


MANUFACTURER DATA 


10490 PARTICLE SAMPLE 


Particle Count (Number } 


0.0 5600.0 1000.0 1500.0 2000.0 2500.0 3000.0 3500.0 4000.0 





0.0 2.0 . 4.0 : 6.0 9.0 10.0 12.0 
Porticle Diometer (Pixels) 


PictncmCn me niulleManutacttirer s Size Distribution 
White Light Reconstruction 


Thus, a modified case is needed. Referring back to the 
Barticle image data in Figure C.2, it is seen that all 
particles smaller Ehiaiileleeeoor Microns gecouwld One be 
reconstructed in the white light case. Figure 4.9 of the 
manufacturer's data has been replotted in Figure C.4 to show 
this threshold somewhere between diameters of 8.89 microns 


pease ill./3 microns. As the particles on the array have 


Go 


discrete sizes, the exact location of this thresho la lie 
no taknown. however, it is known that it lies between these 
two values. The particles which fall below this threshold 
line must be truncated from the manufacturer's data in order 
to accurately predict what can be expected from the IBM/AT 


Sec e ell. 


2500.0  %3000.0 


MANUFACTURER DATA 


10490 PARTICLE SAMPLE 


2000 .0 


THRESHOLD 


Porttcle Count (Number) 
1000.0 1500.0 


$00.0 


0.0 





0.05.0 10.8 15.0 20.0 2.0 %.0 35.0 40.0 45.0 $0.0 $5.0 60.0 65.0 70.0 75.0 80.0 85.0 30.0 & 
Particle Diameter (Micron) 


Figure C.4 Manufacturer ’5S Data Showing 
White Light Size Threshold 
Truncating the manufacturer's data at this dtamemems 
threshold, the numbers shown in the final column Of oi aoe 
C.1 are found. It can be seen that all of the one pixel 


particles are excluded, as are some of the two pixel size 
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Particles. Tiiicmetmiimesyted Smanutacturer s data has been 
plotted and is shown as Figure 4.10. 

The analysis and conversion of the manufacturer's data 
HOLE ne = NoLOgsram reconstruction will now follow the same 
EmomMcmasethose taken im the white light reconstruction. 

Again in this case, an image of the calibrated array was 
Faken at @ 2X magnification and processed. The method of 
filtering and the level of thresholding was recorded so that 
it could be reproduced for the reticle image later. A 
photograph of the hologram array is shown in Figure 4.6. The 
image was then stored and the Fortran program was run_ to 
extract the feature data. EPWestabulartwoutput is shown in 
Pee Ur ere 6. >. 

Uictege ew wrncm.ewidth data as shown in Figure C.5, for all 
of the particles which were present and which were obviously 
not malformed, an average diameter per pixel was calculated. 
The diameters in microns for each of the particles was taken 
meom Ene Manufacturer s data in Figure C.1. The value for 
PeicGhuumGeatUure Sis shown in Figure C€.5, as is the total 
Beverage. (it should be noted that the data for particle / 
was not used in forming the average length per pixel as the 
data had too large a variance to be considered.) 

This average diameter per pixel value was then again 
Msed as the conversion factor to convert the manufacturer's 
Dabticle distributions to pixel diameters instead of micron 


diameters. 


TD 


PART UGE AREA ERE OID IL st Y WIDTH "DIA/ Pies REMARKS 


NO. (PIXEGS) (PIXELS) (CPEXERSe 
al MISS ie 
2 MISSING 
5 MISS INe@ 
4 MISS Bile 
> MISS iNe 
6 Miss Die 
yh 1 ik ii Zon OU NOTE (ae 
8 2 jue i 13.33 
9 3 2 2 i i Jaren 
10 4 ZL 2 17 20 
i | 6 2 5 16 s2oKG 
lez 8 3 4 1325.6 
AS TO 3 3 14.40 
14 L7 4 =) Lies 
1 19 4 5 | iy Ae ps 
16 14 4 5) 13308 
5 Ie ZS 5 6 ier 2 7 
1S 21 5 6 12. 
19 24 6 6 akon F200, 
ZU 42 7 9 1 Org26 
ah 40 7 8 11.204 
Zee ah 8 8 LO. 
ZS 48 8 8 ge arts a 
AVERAGE = 13.08 MICRONS /F baa 
NOLES TGiG) Particle not included in total average 


as it has too to large a variance 


Figure C.5 Calibrated Arrays eatiuwesviane 
Hologram Reconstruction 
Using the conversion factor calculated as shown, a table 
for the hologram reconstruction was made. The table shows 
the width range in pixels verses the width range in microns 


and the total number of particles which fall into thaGeleaee 
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TABLE €.2 CONVERTED MANUFACTURER'S DATA 
HOLOGRAM RECONSTRUCTION 


Diva t ER DIANE LOR TOGRAL NG. OF ~PRUNCATED NO. 


(PIXELS) (MICRONS) PARTICLES OF PARTICLES 
Ce jf) SS) Gal 5749.9 0.0 
iL cay Loi = 2a? 2608.1 5 E15 
es Ge 2 me. 2 988.7 988.7 
= On 2 ame? 3 738.0 738.0 
4 - 5 ye e654 298.0 298.0 
5 - 6 Sgt Sess 90.3 90.3 
6 - 7 moo) = On. 5 16.6 16.6 
7 - 8 91.5 - 104. ee0 iL 


6000.0 










$000.0 


MRNUFACTURER DATA 


10490 PARTICLE SAMPLE 


2000.0 300.0 4000.0 


Particle Count (Number ) 


1000.0 


0.0 


0.0 2.0 4.0 6.0 9.0 10.0 12.0 
Porticle Dicmeter (Pixels) 


Ppcremoro ms Pull! Manufacturer s Size Distribution 


Hologram Reconstruction 
a 


vi 


range. The results are shown in the first three columns oe 
Tape rc. 2. 

A histogram of the total manufacturer's size 
distribution based on pixel diameters has been plotted and 


ismshown in FigcuresG. 6. 


3000.0 







2500.0 


MANUFACTURER DATR 
10490 PARTICLE SAMPLE 


2000.0 


THRESHOLD 


Particle Count (Number) 
1000.0 1500.0 


$00.0 


0.0 


6.05.0 10.0 15.0 20.0 25.0 30.0 35.0 40.0 45.0 50.6 $5.0 60.0 65.0 70.0 75.0 60.0 85.0 90.0 & 
Particles Diameter (Micron) 


Figure C.7 Manufacturer's Data Showing 
Hologram Size Threshold 
Again a modified case is needed to remove all of the 
smallest particles which were not reconstructed in the 
hologram case. Referring back to the particle image data in 
Figure C.5, it is seen that all particles smaller than 24.00 
microns could not be reconstructed. Figure 4.9 of )3Gme 


manufacturer's data has been replotted in Figure C.7 to show 
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this threshold size line somewhere between particles of 
diameters of 21.46 microns and 24.00 microns. The particles 
which fall below this threshold line must be truncated from 
the manufacturer's data in order to accurately predict what 
can be expected from the IBM/AT system. 

Truncating the manufacturer's data at this diameter 
eheeshold , the numbers shown in the fourth column of Table 
C.2 are found. It can be seen that all of the one pixel size 
particles are again excluded, as are some of the two pixel 
size particles. This truncated manufacturer's data has been 


plotted and is shown as Figure 4.12. 


(ee 


AP PEN Di Xap 
CALIBRATED RETICLE EXP ERENT Aleta, 
WHITE LIGHT RECONSTRUCT re, 

The calibrated test reticle was reconstructed wsaing 
white I 2 hte It was then filtered. processed and 
thresholded in the same manner as was the white light 
illuminated calibration test array, shown in Figure 427 Gage 
Fortran particle data extraction program was then run on the 
Converted" .Gata: sfaidbe. The histogram data table from. the 
program is shown in Table D.1 and the tabular owWtp Wc 


shown in Table D.2 


TABLE D.1 HISTOGRAM FEATURE DATA OUCTEUE 


TOTAL NUMBER OF FEATURES = "152 
TOTAL NON-ROUND FEATURES = 7 
BEN sSilZee Ge Lxens) NUMBER OF FEATURE 
FEATURES DENS Tie 
FROM GREATER TO LESS 
THAN THAN 
Or 40 i) 14 0.109) 
se, ite) a2 Ol 2a 
iE EN, oe G a8) Ons 
3.0 4.0 12 O be 
4.0 3.0 10 Oe. 
52.0) 6.0 gles 0:5 ae 
O-.0 Vad 7 0.05 
Hoe) S20 § 0202 
82.0 oF0 1 Os0n 
9.0 16107510. 3 Orm0r 
10.20 1 ee ©, it O20 
ee sa, 2 O20 
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TAD Eee LABDUDARM2ALURE = DATA OUTPUT 


FEATURE NO. AREA X WIDTH Y WIDTH REMARKS 
i 8 3 3 
2 21 4 6 
3 4 2 2 
4 8 3 3 
5 14 4 4 
6 72 10 10 
7 19 6 5 
8 i 4 4 
9 10 3 4 

10 2 1 2 
i 6 3 3 
ie 7 2 4 
LS 36 6 7 
14 4 2 2 
65 44 6 8 
16 10 3 4 
ie7 88 iL 10 
18 24 6 6 
19 iS 4 5 
20 28 6 6 
21 4 2 2 
An 58 7 10 NOT ROUND 
23 4 2 2 
24 7 3 3 
25 18 4 5 
26 8 3 3 
27 2 2 1 
28 12 4 4 
29 BIS 6 7 
30 4 2 2 
34 6 2 3 
32 5 3 3 
33 3 2 2 
34 8 3 3 
35 1 1 1 
36 6 z 3 
37 10 3 4 
38 6 2 3 
39 34 6 7 
40 ie 3 4 
41 7 3 3 
42 7 3 4 
43 6 3 3 
4d 4 2 2 


ol 


FEATURE N@@ AREA X WIDTH Y WEDTH REMARKS 


45 8 2 5 NOT ROUND 
bo 19 mn 6 : 
47 29 6 6 
48 29 5 i 
49 7 3 3 
50 7 3 3 
51 i i fl 
52 32 8 6 
53 30 7 6 
54 r 2 2 
55 26 5 6 
56 6 2 3 
57 ath 5 5 
58 64 el 7 NOT ROUND 
59 8 3 4 
60 Di 5 5 
61 d 2 2 
62 23 5 5 
63 74 9 10 
64 i 2 3 
65 2 1 2 
66 29 6 7 
67 14 mA 5 
68 16 d 5 
69 7 3 3 
70 7 4 5 
a 9 3 d 
72 14 mA 4 
73 34 6 7 
74 d 2 2 
75 i 2 2 
76 1 1 1 
77 34 7 7 
78 14 4 5 
79 i 2 2 
80 1 1 1 
81 2 2 1 
82 a 5 6 
83 63 iB 9 
84 ais 6 7 
85 12 3 mA 
8 6 ha 3 4 
87 5 3 3 
88 a5 6 7 
89 6 2 3 
90 2 1 2 
91 35 6 7 
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FEATURE NO. AREA X WIDTH TewlDta REMARKS 


92 116 ie 14 
93 19 4 6 
94 3 2 2 
95 ial 3 5 
96 18 5 5 
97 2 2 1 
98 y 3 3 
99 7 3 3 
100 2 fl 2 
ion iL 7 4 5 
102 3 2 2 
103 10 4 3 
104 6 2 3 
105 54 7 9 
106 56 8 8 
107 5 3 2 
108 33 6 i 
109 17 4 5 
110 34 6 7 
ial 6 2 3 
ii? fl 1 1 
113 48 7 8 
a 20 4 6 
Las 7 3 3 
116 7 3 3 
Lily 95 10 2 
iL al fe : 58 10 7 NOT ROUND 
119 9 3 4 
120 1 1 1 
121 2 1 2 
fh22 10 3 4 
123 8 3 3 
124 8 3 3 
125 37 6 7 
126 1 1 1 
27 19 5 5 
18 46 7 8 
129 12 3 4 
130 19 5 5 
131 4 2 2 
fe3)2 2 2 1 
133 iil 3 5 
134 54 8 9 
LSS 33 6 7 
136 4 2 2 
137 rae 8 18 NOT ROUND 
138 39 7 7 


33 


BEATURE NO. AREA X Wea Tawa Died REMARKS 





SS) 7 3 3 
IL 8) 2 al Z 
iL ak 6 Z 3 
ea 48 yi 9 
43 124 15 Dol NOT ROUND 
144 19 4 6 
eye 26 5 6 
146 64 120 8 
147 4 2 2 
Tas 10 4 5 
149 8 > Z NOT ROUND 
15:0 Al 1 i 
jist 2 Z i 
Note.c.) The data from these tables is shown in graphical 


Porn un Figure 4.271) 
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PEPENDIX E 
COP oR weE Dek ETECLE EXPERIMENTAL DATA 
DASER HOLOGRAM RECONSTRUCTION 
The calibrated test reticle was reconstructed using 
the laser hologram. It was then filtered, processed and 
thresholded an the same manner as was the laser 
Pecomnstructed calibration hologram test array, shown in 
Peeure 4.6 . The Fortran particle data extraction program 
was then run on the converted data file. The histogram data 
table from the program is shown in Table E.1 and the tabular 


Smtput is shown in Table E.2 . 


Tabu wie HESTOGCRAM FEATURE DATA OUTPUT 


mera NUMBER OF FEATURES = 228 
TOTAL NON-ROUND FEATURES = 11 
Biles i ZESCEILXELS ) NUMBER OF FEATURE 
FEATURES DENS 
meoM GREATER To SES > 
THAN THAN 
0.0 sre 47 OG fagll 
iO 20 60 O26 
220) S00) Bue OO. 25 
B70 4.0 30 Om 
4.0 20) eZ ORO: 5 
a0 OU 6 ORCS 
G20 Ua”) 2 Oe0 1 
JO 6700 Z Ok Ola 
oie) 3) 0) 0 02010 
2 <0 TO) ee 0 05:00 
io0 Ia © 0 O00 
ge 8 20) 0 O00 


oe 


TABLE E.2 TABULAR FEATURE SDA Tc Tere 


FEATURE NO. AREA X WIDTH Y WLPTH REMARKS 


NOT ROUND 


CONDOM” fF WN FP 
a | 


NMNNPP PRP PRP PRP PB 
NF OOMNI AUP WNrR OO 
k 


NOT ROUND 
NOT ROUND 


BR bo fh he 
Om f& W 
nm FPR 


EP PW WW WW WW WD WWD PH 
NP ODAN AUF WHF OO OO 
ho _ = 


i 
i> 
MP PHN RP NPNHKP ER FWEN UN WR UF ADANBHPWNHWNW KP RPNWHRN FN UDWHDND F&F 


MmMrPrP FW EWM WNW AW UkP FM WR AIP WOrR WW EH FWRPR NWR DY FH WMO WH DS ™“ 


fo 
~~ 
WRPRPAFUWWUN OWMINUNAFNRP PRP UWPNAI DOK AOArPPRPPWOPR 


> 
> 
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eA Lire NO. AREA X WIDTH Y WIDTH REMARKS 


i al 


45 
46 i 
ae 
48 
49 
50 
el 
az 
aE. 
54 
55 
56 
3) 
Bae. 
3 Ht 
60 
61 
62 
63 
64 
65 
66 
67 
68 2 
69 
70 
al 
72 
ies 
74 
i 
76 
Tied 
Tae: 
79 
80 
81 
82 
2 
84 
85 
86 
a7 
88 
39 
20 
sal 


NOT ROUND 


ead 


—_— 


co 


js 
Wo FP OrRrN UWA NM FF ORF WON DAF Fr WWRrR OW OND FW YY 


NOT ROUND 


lo 
ho ON 
ray 


a — 
OUU AAO 
WNHRPWREFPWNHE RP RPWWUDONNP EF RPRPAYINNNWERWWWUNWEPNPFP UP RP FH PPP RP 


PmrWWH WR WW UND WWI DN WP WRN UWP NM NW FW EF WWPRPNP FFP WN WH LF UW bd 


OF rPn~ WwW BH 4 DB ~™ 


eal 


OW 


PEALURE NOe AREA X WIDTH Y WIDTH REMARKS 





ge : 2 
73 
94 
og 
70 
a7 
oe: 
79 
LOO 
On 
OZ 
LOS 
Oz 
POS 
106 
LO7 
108 
109 
110 
ileal 
1 EY 
pap les, 
114 
ia e 
116 
ey 
Tae 
ile, 
1720 
Eval 
22 
23 
124 
2 
WZ 
ili 
aye 
PZ9 
130 
Lt Z 
AS 2 
33 
134 
S|) 
6 
ies, 
aL Site: 


NOT ROUND 


F* fo 


0o 
Mor RP WwW WOW OUWOWR fF UOUWOnNRPN PE WOR WW 


oe 


—_ 
Ne) 


NOT ROUND 
NOT ROUND 


ao eo 
UWrR~ Fr Oo 


ON 
FPtoW RPM WH WWW fF WHR DY FF FWWWE DN WNP WKH WH UF HY WNP HNN FWWOWHNY MH WEF KH “J 


mM fF FP RPM WWW WWNHR WEF FWE ON WF WR WHE HM UUND ONR NY WH WWW EWN HK WENDY F 


Pro OO RP KH UND “4 On WH 
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FEATURE NO. AREA X WIDTH Peewee DTH REMARKS 


spe, 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
0 
eo 
Ie 
Ia 
154 
ipo, 
156 
Lo? 
ae 
ils s, 
ioe 
Isa 
Gi 
163 
164 
Gn 
166 
16 7 
Ihisys 
iG 
Ly) 
eed: 
ie 2 
73 
174 
i as 
NAF? 
177 
IL iis 
iy 9 
mone) 
31 
eS 2 
33 
184 
ik reyes, 


me 


ho 


ho bo 


NOT ROUND 
NOT ROUND 


a 


ho 


re 
COW WW FMAM BD WR WWNH NHR NMNM NH MH DH HW FN RrP RP RP FF WWNHN DOP OW WW FR WR PRP 


DWWUNYUNP FPN EN FPN NYNNHNH NN FSF WEP PHY NWWWNHNAUKPDWNHNM WP UNEP 


eo = 
IMO PfwFr FwmoworOnwaEFarRrFPEF WE WEEDON PRM AYNOnA FWE HR WWE DOP FHP LF 
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PEATURE NO AREA X WIDTH Xt WEDTH REMARKS 





66 J) ik Z 
Jie) 7 7 3 i) 
isis: > 5 2 
189 8 5 s) 
1990 Ay 4 5) 
1 oat 9 3 4 
192 Ais. 3 5 
je ge 18 3) 2B) 
194 4 2 2 
195 3 2 2 
196 10 3 4 
197 10 $ 4 
198 LES) 3 4 
199 Zo 4 6 
200 28 10 3 NOT ROUND 
2 Ol 7 2 4 
ZO02 9 S 4 
OS eh 4 4 
204 a d. ak 
205 dies 5 4 
706 4 2 2 
BOF ie i i 
20S 10 eS 4 
209 1 il 1 
ZAG 2 1 2 
pag | P22 4 4 
2 ee 1 Al i 
yg BS 14 4 4 
ZAR 3 Z $ 
ZS 2 1 Zz 
Z16 6 3 3 
Bealed 4 2 2 
21.8 i 4 6 
219 4 3 a 
Za aE 1 il 
Z21 5 Z S 
(a hes 1 4 i 
223 3 iz Z 
224 4 J) Z 
225 es 6 4 
Dass: ay, 5 4 
227 6 5 3 
PPK a al ak 1 
Nowe.) The data from these tables is shown in graphical 


form in Pigure,4.44 
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